Entropy production at freeze-out from dissipative fluids 
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Abstract: Entropy production due to shear viscosity during the continuous freeze-out of a lon- 
gitudinally expanding dissipative fluid is addressed. Assuming the validity of the fluid dynamical 
description during the continuous removal of interacting matter we estimated a small entropy pro- 
duction as function of the freeze-out duration and the ratio of dissipative to non-dissipative quantities 
in case of a relativistic massless pion fluid. 
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I. INTRODUCTION 

In the fluid dynamical modeling of high energy heavy- 
ion reactions the initial large gradients dilute and cool 
the matter until the mean free path of particles becomes 
comparable to the size of the interacting system. By then 
the particles do not interact strongly enough to maintain 
local equilibrium and the interactions between particles 
cease leaving the momentum of particles frozen-out [l[ . 

To properly describe the final break-up we need to 
model the continuous transition from an interacting fluid 
to a noninteracting fluid (gas) of particles. Here we will 
use the standard first order theory of relativistic fluid 
dynamics, i.e., the relativistic Navier-Stokes equations, 
and calculate the entropy production during the time 
the matter is continuously removed from the interact- 
ing fluid. In reality the matter equilibrates in finite time, 
where the local relaxation time, r t h, increases in time as 
the system dilutes and the interacting matter turns grad- 
ually into a weakly interacting dissipative fluid before the 
momentum of particles freeze-out 0, Q • Assuming that 
the fluid is near local equilibrium during freeze-out we 
can estimate the additionally produced entropy from the 
equations of dissipative fluid dynamics. Strictly speak- 
ing this transition does not necessarily lead to kinetic 
freeze-out, since the mean free path of the matter is still 
finite and the particles may interact further. The drain 
from the fluid should be followed by a kinetic description 
where this transition can be treated beyond small devia- 
tions from local equilibrium, but we will not discuss the 
kinetic modeling here. 



II. FREEZE-OUT FROM DISSIPATIVE FLUIDS 

A single-component fluid is characterized by a con- 
served charge four-current, _/V M , the energy- momentum 
tensor, T^, and the entropy four-current, S^. The dis- 
sipative currents are expressed explicitly by projecting 
the above quantities perpendicular to the flow of matter 
[U, [H] • Hence the conserved charge four-current is decom- 
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posed using the transverse projector, A 
where is the fluid dynamical flow velocity and g 



diag(l, —1,-1,-1) the metric of flat space-time, thus 

= nu» + V , (1) 

where n is the local rest frame invariant charge or con- 
served particle density, and V 11 = A^ N„ is the flow of 
charge perpendicular to . 

In case of a relativistic massless pion fluid we can ne- 
glect the heat conductivity and the bulk viscosity below 
the critical temperature since both quantities are very 
small or vanishing compared to the shear viscosity [1,0]. 
Thus, the energy-momentum tensor is decomposed as: 



(2) 



where e and p are the local rest frame energy density 
and isotropic pressure and the symmetric and traceless 
part of the energy-momentum tensor is the stress tensor, 



±A^A n 



"| (A^ + A^- f ^- 

These conserved quantities satisfy the following conti- 
nuity equations written in differential form, 



d^N* = 0, 
d u T^ = 0, 



(3) 
(4) 



where, <9 M = d/dx^, is the space-time four-divergence. 
These equations have to satisfy, in addition, the condition 
of positive entropy production required by the second law 
of thermodynamics, thus the entro py p roduction with 
shear viscosity becomes % i, ES El El El : 
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> 0. 



(5) 



Here the shear stress tensor is, n^ 1 ' — 2rja^ u , where r\ 
is the shear viscosity coefficient and the stress tensor is, 
a*"* = \ {Wu u + V^)-±A^V a m a , and V = A»"d v . 
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In the perfect fluid limit, (jth 
production is zero [3, EB| : 



0, 7? — ► 0), the entropy 



^ er/ = -^(n<) = 0, 



(6) 



where T is the temperature, and fi is the chemical po- 
tential. The above equation is satisfied if the number of 
charge(s) is conserved or for a vanishing chemical poten- 
tial. Aiming to study the freeze-out of massless pions the 
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later condition is favored since fi = [lq . 

To describe the freeze-out of matter, one has to sepa- 
rate the fluid into two components [lj]: interacting and 
frozen-out, (denoted by lower indexes i and /). The in- 
teracting component continuously emits particles and its 
evolution is described via the equations of fluid dynam- 
ics. The emitted particles leave the fluid and part of it 
may interact further described by Boltzmann transport 
equation [ll| , until the particles will no longer collide 
but free stream from their kinetic freeze-out points. 

If the particle number is conserved during freeze-out, 
d^N^ — 0, we can separate the system into interacting 
and frozen- out particles, = JVf + Nf, such that, 



(7) 



Similarly, but independently (even for non-conserved 
charges), the total energy and momentum is conserved, 
d^T^ = 0, thus the change in the energy-momentum 
tensor at freeze-out leads to: 



-r. 



-J" 



(8) 



where, J u — (J , J), is the drain term which takes away 
energy and momentum from the interacting component 
and transfers it to the frozen-out component. Generally 
there is a non- vanishing momentum transfer, i.e., there 
is non-zero mechanical work done on the interacting fluid 
at freeze-out. The equation for the entropy, constrained 
by the second law of thermodynamics is: 



(9) 



where the entropy production in the frozen-out matter is 
assumed to be zero. 



1. Particle number conservation, energy conservation and 
entropy production 

In what follows, we will work in the light-cone coordi- 
nates, making the coordinate transformation from (t, z) 
to (t, rj) , where r = y/W 



2 is the proper time and 
r) = |log((i + z)/(t — z)) is the space-time rapidity. 
Here we present the equations of dissipative fluid dy- 
namics with a drain term using the 1+1D boost invariant 
Bjorken model [2(J, where m m = (cosh 77, 0, 0, sinh^) is the 
flow velocity and 9 = = 1/r, is the expansion rate 

independent on the equation of state. 

Assuming that the particle number is conserved leads 
to the following equation for the interacting density: 



dn t (r) = [-dNi(r) - ni(r)dVi(T)} /V^t) , 



(10) 



where the volume of the interacting matter Vi (r) ~ tAt 
increases such as, c?V^(r)/V^(r) = oIt/t, and the trans- 
verse area of the system is At ■ The volume of the frozen- 
out particles is unknown from fluid dynamics but it is of 
no relevance whatsoever, hence in the following we will 



not discuss the frozen-out component. 

Furthermore, for simplicity we assume that the change 
in the number of interacting particles is due to parti- 
cle freeze-out only, governed by a boost invariant proper 
time dependent probability, P esc (r)dT , such that: 



dN l (T)=-P esc {T)N l ( T )dT. 



(11) 



where dNf = —dNi and P esc (r) is the escape rate. The 
properties of matter change considerably in a small space- 
time domain, characterized by a time-scale, t'(t), which 
governs the process. Here, we recall the assumption for 
the escape rate [U [22[, which characterizes the freeze- 
out of matter in finite time, L, the proper time interval 
between two parallel hypersurfaces, where the continuous 
freeze-out of matter starts at tq and ends at tq + L, thus 



Pesc{T) = 



t'(t) t'(t ) L + t -t ' 



(12) 



where the parameter, t'(tq), denotes the initial charac- 
teristic freeze-out time, see Sect. IIII1 For example at the 
midrapidity section of the hypersurface, the probability 
of escape depends on the invariant distance x M d 3 cr AI ~ r 
in the direction of the freeze-out normal vector and on 
the local properties of matter, such as the proper den- 
sity and invariant cross section, i.e., the collision rate, 
r ~ na. Therefore, t'(t) ~ T, is inversely proportional 
to the relaxation time [H, [22j ■ Fluid dynamics re- 
quires that r ^> 9, while the kinetic freeze-out condition 
is T ~ 9 [H, [24], HBj] . Now, using the above equations, the 
change in the interacting particle number density due to 
simultaneous Bjorken expansion and freeze-out becomes: 



dr 



= -P esc (T)rii(T) 



ni(r) 



(13) 



Let us recall the first law of thermodynamics neglecting 
the baryon number and strangeness conservation, that is, 
dEi = TidSi — PidVi, where PidVi denotes the work done 
by the expansion. Here, TidSi = AQi+TidSf lss denotes 
the change in the total entropy at constant temperature 
and TidSi > AQi expresses the second law of thermo- 
dynamics, where the entropy increase, TidSf lss , is due 
to dissipation (shear viscosity). Thus, for an adiabatic 
expansion, AQi = 0, the change in the total energy is: 



dEk(r) 



[ Pl (T)-^( T )}dV(T) 



(14) 



where Pi(r) is the isotropic pressure and $i(r) is the 
stress which slows down the expansion and decreases the 
"useful" work done by the pressure. Furthermore, the 
pressure is independent on the rapidity, dpi(r)/drj = 0. 

Now, let us consider simultaneously the adiabatic ex- 
pansion and the freeze-out of the system. The decrease 
in the total energy of the system is due the energy taken 
away by the escaping interacting particles and the work 
done on interacting component by the frozen-out compo- 
nent during the particle escape from the freeze-out sur- 
face [26|] . With a similar ansatz as in eq. (fTTj) we obtain 



3 



an equation for the evolution of the energy density: 



dr 



-[ei{T)+ Pi (T)]P esc (r) (15) 
ei(r) +Pi(r) $,-(t) 



Furthermore, using the first law of thermodynamics we 
get equation for the evolution of the entropy density: 



dr 



-Pesc{T)s l {T) 



Si(r) $i(r) 



Ti(r)r 



(16) 



Here the shear is defined to be positive, <f> = it"" — n" 
In the standard first order theory, $ = <S>i st = (4?y)/(3T), 
where r; is the shear viscosity coefficient. In second order 
theory of relativistic fluid dynamics [27| , the initial value 
for the shear and its evolution must be obtained indepen- 
dently. Introducing the Reynolds number: R = (e+p)/$, 
which represents the ration of non-dissipative to dissi- 
pative quantities (28l . [29j . then we can choose the ini- 
tial value for the stress such that dR(r )/dT = [30| . 
Even in second order theory this physically motivated 
initial value for the stress brings the behavior of the sys- 
tem very close to the standard first order theory during 
the evolution, but in any case the <E> relaxes to the first 
order-theory, given by the relaxation time for the shear. 
Therefore in the following we use the first order theory to 
describe the late stage of the fluid dynamical evolution. 



III. RESULTS AND DISCUSSIONS 

We assume that we deal with a massless pion fluid at 
freeze-out with an EOS of an ideal fluid, e(r) = 3p(r), 
where the energy density is, e = aT 4 , a = g Tr 7r 2 /30, g n = 
3 is the degeneracy of pions, while the entropy density 
is s = 4aT 3 /3. The particle density is n = 6T 3 , where 
b = g 1T C(3)/TT 2 . Here we recall the result of Gavin 0] for 
the shear viscosity coefficient, (see also [|[ and references 
therein), in case of a massless pion fluid, calculated in the 
relaxation time approach from the relativistic Boltzmann 
transport equation, 



nvn = 15 e(r) 



(17) 



The EOS together with eq. (H]) and eq. fT7). gives the 
following differential equation for the temperature, 



dT ^ X p ( W f 1 



/ tT th {T ) \ dr 
45r 



(18) 



leading to the following solution with the escape rate eq. 
(fT2)) for fixed relaxation time, i.e., T{h(r) = t^to), 



T(t) = T 



T -T 



3t'(to) 



where T = T(tq). The result for a perfect fluid can 
be obtained in the limit when, r t ^(r) — > 0. Here we 
also take into account, with a simple ansatz, that dur- 
ing freeze-out the relaxation time increases as the sys- 
tem dilutes and becomes sparser because of the expan- 
sion and the removal of interacting matter. Thus even 
if the freeze-out starts from a perfect fluid, the relax- 
ation time increases during the finite time, L, such that: 
T~th{T~) — > Tth(To)(r — tq)/L, reaching the dissipative fluid 
limit, T th (T ), at r = r + L. 

For the purpose of illustration, and to obtain numerical 
results, we assumed the following initial values for our pa- 
rameters. The continuous escape of particles lasts, L = 4 
fm/c, measured from the start of the freeze-out, tq = 5 
fm/c. We assume that the initial thermalization time 
and freeze-out characteristic time is r t ^(ro) ~ t'(to) ~ 1 
fm/c. This is consistent with the average mean free time 
between collisions estimated from two-particle correla- 
tions [3l|. The initial temperature is To = 160 MeV. 

The evolution of the temperature is shown on Fig. [T] 
Here T\ and T-i show the decrease of the temperature of 
the interacting matter during freeze-out from a perfect 
fluid and dissipative fluid respectively. The latter cor- 
responds to eq. (fT9|) . while T\ corresponds to the same 
equation in the perfect fluid limit. The temperature de- 
creases to zero in both cases signalling clearly the com- 
pleted removal of the interacting matter. 




(19) 



x - x Q [fm/c] 



FIG. 1: The temperature of the interacting matter as a func- 
tion of the freeze-out duration t — To- Ti is for a perfect fluid 
and T2 is for a dissipative fluid. 

We have also employed the ansatz for an increasing re- 
laxation time, Tthij) — Tth{To)(r ~ tq)/L, and checked 
our result against this approach; the temperature curve 
is essentially the same as earlier. In all cases the viscos- 
ity is small and does not affect the fast decrease of the 
temperature due to the freeze-out of matter, but in case 
of large viscosity the temperature decrease is slower than 
in the perfect fluid limit. 

Fig. [2] shows the ratio of shear viscosity to entropy 
density and the inverse of the Reynolds number, using a 
fixed and a dynamically increasing relaxation time. From 
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eq. (fl~7|) the ratio of shear viscosity to entropy density is: 
i](t)/s(t) = T(r)rt/j(r)/5. Here r?i/si shows the ratio of 
shear viscosity to entropy density for a fixed relaxation 
time. Initially this value is close the conjectured minimal 
bound, 77/s > 1/(4tt) l32l ■ (see also 0,111], and in fluid 
dynamical simulations [13l l34j|). In this case, f?i/si is 
decreasing with the temperature. The next, 772/ 'S2, ratio 
employs the increasing relaxation time. Since the mean 
free time between the collision increased so did the abil- 
ity to carry momentum, i.e., viscosity, and the ratio in- 
creases [33| . Since we remove energy (momentum) from 
the system, the remaining cold matter has less energy 
(momentum) available to transfer causing the decrease 
of the ratio. 



0.2 




x - x Q [fm/c] 



account through the increasing thermalization time, the 
total entropy produced at the end is as much as during 
the evolution of a dissipative fluid undisturbed by the re- 
move of interacting matter, since the above ratio is the 
same with or without the freeze-out term. 
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x - x Q [fm/c] 



FIG. 3: The entropy production as a function of the freeze-out 
duration r — r . Here Si and S2 corresponds to the entropy 
production at the freeze-out from a dissipative fluid with con- 
stant and increasing relaxation time respectively. 



IV. CONCLUSIONS 



FIG. 2: The ratio of shear viscosity to entropy density, and 
the inverse of the Reynolds number as a function of the freeze- 
out duration r — r ; see the text for more explanation. 

The inverse of the Reynolds number is: i? _1 (r) = 
4rifc(r)/(15r). Here, R^ 1 shows the evolution of the 
ratio of dissipative to non-dissipative quantities with a 
constant relaxation time. R^ 1 shows the transition from 
a perfect fluid to a dissipative fluid during freeze-out, 
where the viscosity of the remaining interacting compo- 
nent is increasing with time. 

The next figure, Fig. [3l shows the entropy production 
as a function of the freeze-out duration, for a fixed re- 
laxation time 5*1 and for the increasing relaxation time 
S2 respectively. We can explicitly estimate the entropy 
production comparing the evolution of the entropy den- 
sity, s = 4aT 3 /3, of a viscous fluid to a perfect fluid and 
calculate the ratio of entropy densities, s dlss / s P er f : thus 
for a fixed relaxation time, 

^l=expf^fl-i)l. (20) 

The entropy production is a function of the freeze-out 
duration and of the ratio of dissipative to non-dissipative 
quantities, R , but it is independent on the continuous 
removal of interacting matter. Thus during the grad- 
ual transition from perfect to dissipative fluid, taken into 



We have incorporated a simple freeze-out criterion into 
the fluid dynamical model, which continuously removes 
the conserved charge(s), energy, momentum and entropy 
from the system. Since, entropy production is related 
to the speed of equilibration while the freeze-out drives 
the system out of equilibrium, it is also important to 
include and treat the deviation from equilibrium of the 
remaining interacting matter as this produces additional 
entropy which should be added when estima ting the to- 
tal entropy produced in heavy- ion collisions [30T[35l. |36| . 
As long as these two processes are comparable and the 
deviation from equilibrium is small, the entropy produc- 
tion is given by the standard viscous correction to the 
equations of fluid dynamics. 

The result of our modest and illustrative analysis 
shows that the produced entropy is small even for long 
freeze-out times. Compared to the entropy produced in 
a short time, i.e., ~ 1 fm/c in the initial stage, during 
the same amount of time at freeze-out this corresponds 
up to two orders of magnitude less entropy production. 
Here one could use different relaxation times, a different 
EOS and take into account higher order corrections in the 
entropy four-current to (de-)increase the entropy produc- 
tion, but still in case of a massless pion gas this would 
not be large. Finally, taking into account the transverse 
expansion the entropy production due to shear viscos- 
ity would slow down and by the time the longitudinal 
and transverse size of the system become comparable, 
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the shear viscosity vanishes [6]. However, even in this 
case substantial entropy production can be expected in 
a hadron gas where the transport coefficients are differ- 
ent while the bulk viscosity and heat conductivity are no 
longer negligible. We believe that such calculation has 
much to clarify on the question of entropy production in 
high-energy heavy-ion collisions. 
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